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Abstract 



We perform the NLO QCD analysis of the world data on inclusive deep 
inelastic scattering cross sections of charged leptons off the proton and the 
deuterium targets. The parton distributions, the value of strong coupling 
constant as, and the twist 4 contributions to the structure functions F2 
and Fl are extracted with the complete account for the correlations of 
data points due to the systematic errors. Sensitivity of the value and 
the high twist contribution to the procedures of accounting for the sys- 
tematic errors is studied. The impact of theoretical uncertainties on the 
value of as and on the parton distributions is analysed. The obtained 
value of strong coupling constant with the account of these uncertainties is 
as(Mz) = 0.1165 ± 0.0017(stat syst) ±'^S (theor). The uncertainties of 
parton-parton luminosities for the FNAL and LHC colliders are estimated. 



PACS number(s): 13.60Hb,12.38Bx,06.20.Jr 
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I. INTRODUCTION 



Experiments on deep inelastic scattering (DIS) of leptons off nucleons is a unique source 
of information about strong interaction. These experiments were initiated at SLAG linac 
and later were continued at different accelerators using the fixed targets and the colliding 
electron-proton beams. The data for proton and deuterium targets, given in Refs. are 
especially valuable, since no heavy-nucleus corrections are needed for their interpretation. 
Those data combined with the results from HERA electron-proton collider allows one 
to determine the parton distribution functions (PDFs) and are widely used for this purpose. 
In particular, global fits of PDFs, which are regularly updated by collaborations MRST 
and CTEQ [||, rely heavily on the DIS data. It is often mentioned, that the MRST 
and CTEQ PDFs lack information on uncertainties, that does not allow one to estimate 
the uncertainties on the cross sections, which are calculated using those PDFs. Most often 
these uncertainties are estimated as a spread of results, obtained using different PDFs sets. 
Meanwhile, it is evident, that if different PDFs are based on the same theoretical model 
fitted to a similar data sets, this spread mainly reflects uncertainties of calculations, rather, 
than real uncertainties arising from statistical and systematic errors on the data used for the 
extraction of PDFs. Besides, those collaborations combine statistical and systematic errors 
in quadrature, i.e. do not account for the correlation of the latter. Since systematic errors 
dominate over statistical ones for many DIS experiments, they govern total experimental 
errors on the PDFs parameters fitted to the data and ignorance of their correlations may 
lead to the distortion of the parameters errors and to the bias of their central values. 

Statistical and systematic errors are combined in quadrature in part by historical reasons. 
The other reason is that, contrary to the case of statistical errors, existing approaches to 
the account of systematic errors are not so straightforward and encounter with technical 
difficulties generated by correlations between measurements, which become more significant 
when the systematic errors rise, as compared with the statistical ones. Nevertheless, as it 
was shown in Ref. on the example of combined analysis of DIS data from Refs. [0,^,^,^ 
with the complete account of correlations due to systematic errors, these difficulties can be 
overcomed using in the fit an estimator based on the covariance matrix. The results of the 
combined analysis of data from Refs. which attempted to account for correlations of 

systematic errors, was later given in Ref. |T^, but due to the large number of independent 



sources of the systematic errors, they were combined with the statistical errors partially. 
Regardless of the expressed confidence that this procedure should have minimal impact on 
the results, this point is not ultimately clarified and it is evident that errors on the obtained 
PDFs may be distorted. 

In this paper we describe the results of the combined analysis of the world data on 
the charged leptons DIS off the proton and deuterium targets given in Refs. In 
comparison with our previous fit of Ref. ^ in the present analysis we use data with lower 
values of transferred momentum Q. Besides, the data from the SLAC experiments and the 
experiment FNAL-E-665 are added. As well as in Ref. we extract from the data the 
nucleon PDFs and the value of strong coupling constant ctg. In addition, wealth of data at 
low Q allows one to determine the high twist (HT) contributions to the structure functions 
F2 and Fl as well. Analysis is performed in the NLO QCD approximation with the complete 
account of correlations due to systematic errors within approach described in Ref. |jlTI . 
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II. THE DIS PHENOMENOLOGY 



It is well known that the DIS cross section of charged leptons off nucleons can be expressed 
in terms of structure functions -F2,3,lQ- example, at 4-momentum transfers Q less, than 
the Z-boson mass the charged leptons cross section reads 



dxdy 



(2.1) 



where s is the s.c.m. energy, m\ is the lepton mass, y is the ratio of the energy lost by lepton 
to the initial lepton energy, x is the Bjorken scaling variable, M is the nucleon mass, a is the 
electro- magnetic coupling constant. The structure functions F2^^ depend on the variables x 



and Q. Within the operator product expansion |]13| the structure functions are given by the 
sum of contributions coming from operators of different twists. For the unpolarized lepton 
scattering the even twists larger or equal to two contribute only. Thus with the account of 
the twist-4 contribution 

F2,l(x, Q) = F™^(x, Q) + H,,^{x)^^^, (2.2) 
where p^^'^^'^ gives the leading twist (LT) with the account of target mass corrections, as 



calculated in Ref. [Ol: 



^2 (^'<5) = -57i 72 + 6-^7J^2, (2.3) 

{.x,Q) = [x,Q) + -2 + — — (6-2r)J2, (2.4) 

STMC ^ 



where 



5tmc 



_ 2x _ AM^x^ 

CtMC — TT- T — 1 



1 + g2 ' 

and are the structure functions of twist 2. Such approach allows us to separate pure 
kinematical corrections, so that the functions H2;l{x) correspond to "genuine" or "dynam- 
ical" contribution of the twist 4 operators. Note, that the parametrization ( p.2|) implies, 
that the anomalous dimensions of the twist 4 operators are equal to zero, that is invalid in 
general case. Moreover, there are attempts to estimate these anomalous dimensions from the 



^The comprehensive analysis of lepton-nucleon scattering amplitudes, including notations used 
throughout our paper is given, e.g., in review 
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account of the correlations between partons (see Ref. |]T5[)- Meanwhile, in view of limited 
precision of the data, approximation (p.2|) is rather good (see also discussion in Ref. [1T6|). 

The leading twist structure functions can be expressed in factorized form as the Mellin 
convolution of PDFs q with the coefficient functions C: 



1 dz 



— Ql [z,as{nR),Q/lJ,F] qi{x/z,iJ,F) 



(2.5) 



where index i marks the partons species and ctg is running strong coupling constant, 
dependence of PDFs on Q is described by the DGLAP evolution equations |]T7 



Q ^^gg^^ =1]^ ^Pij [2;,a,(/iR),Q//iR]gj(a;/z,Q), 



The 



(2.6) 



and the PDFs evolution is governed by the splitting functions Pij, which in turn depend 
on ctg. The quantities /xp and /zr in Eqs.(|275|) and ( p.6| ) give the factorization and renor- 
malization scales respectively. In the MS renormalization-factorization scheme, used in our 
analysis, these scales are chosen equal to the value of Q usually. The splitting and coefficient 
functions can be calculated in perturbative QCD as series in as- For the coefficient functions 
these series are completely calculated up to the next-to- next-to-leading order (NNLO) [p!8| ; 
for the splitting functions the next-to-leading order (NLO) corrections are known, while for 
the NNLO corrections a limited set of the Mellin moments as well as some asymptotes. 



are available only (see references in |20|). Nevertheless, there are attempts to analyse the 
DIS data in the NNLO QCD approximation with the consideration of the available moments 
only [pl|-|23[], or modelling splitting functions |]2^. Our analysis is performed in the NLO 
QCD approximation with the use of the splitting and coefficient functions in x-space as they 
are given in Ref. . 

The dependence of as on Q is given by the renormalization group equation, which in the 
NLO QCD approximation reac 



1 



1 



as(Q) as(Mz) 



2% \Mz' 



f3 + l/as{Q) 
/3+l/as(Mz) 



(2.7) 



where (3 
102 



/3i 



47r/3o' ^'^ regular coefficients of /5-function: /3o = 11 — (2/3)nf, Pi = 

(38/3)nf, rif is the number of active fermions, which depends on Q. In our analysis 



nf was chosen equal to 3 for Q < nic, 4 for rric < Q < rrih, and 5 for mb < Q < 'n^t, where 



m, 



mb,mt are masses of the c-, b- and t-quarks correspondingly, and when rif changes, 
the continuity of as{Q) is kept (see Ref. |2^ for argumentation of this approach). The 
choice of the quark mass value as the threshold for nf switching is optional. E.g., in the 
analysis of heavy quark contribution to the DIS sum rules, given in Ref. ||28|, this threshold 
is chosen equal to 6.5mc,b,t- Unfortunately any choice cannot be completely justified, while 



^Analogous equations given in Refs. |^J2^ contain misprints, meanwhile, the calculations were 
performed using the correct Eqn.(|2. 
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the dependence of results on the variation of threshold, say in the interval from m^b.t to 
6.5mc,b,t; generates one of the sources of theoretical uncertainties inherent to this analysis. 
Since the value of ctg depends on the threshold position logarithmically, for estimation of this 
uncertainty we shifted this threshold value to the logarithmic centre of this interval, i.e. from 
mc,b,t to \/6^mc,b,t- Very often, approximate solutions of Eg. (|2.7|), based on the expansions 
of as in inverse powers of ln(Q) are used in calculation. Inaccuracy of these expansions for 
evolution of as from O(GeV) to Mz may be as large as 0.001 which is comparable with 



the experimental uncertainties of the as(Mz) value extracted from the data. In order to 
escape these uncertainties we use in the analysis the exact numerical solution of Eg. (^.71) 
instead. 

Since we use the truncated perturbative series, the results depend on the factorization 
scale /ip and the renormalization scale /ir. These dependences cause additional theoretical 
uncertainties of the analysis The accurate estimate of these uncertainties is difficult, because 
the possible range of the scales variation is undefined and besides, one is to change factor- 
ization scheme as well. In our analysis we estimate only the theoretical uncertainty due to 
the choice of /iR in the evolution eguations ( |2.6| ) using the approach described in Ref. pO 



In accordance with this approach the renormalization scale //r is chosen egual to k^Q and 
the NLO evolution eguations are modified in the following way 

= E /; [p!P(.) + ^ [4-)(.) + AP<»'(.) } ,,(./.. Q). 

where p(°) and p(^) are respectively the LO and NLO coefficients of the splitting functions 
series. The change of results under variation of /cr from 1/2 to 2 gives an estimate of the 
error due to renormlization scale uncertainty. Evidently that, by definition, this uncertainty 
is connected with the impact of unaccounted terms of the perturbative series. 

In order to obtain the PDFs from evolution eguations, one is to supply a boundary 
conditions at some starting value Qq. The x-dependence of PDFs cannot be calculated from 
the modern strong interaction theory, it is determined from the comparison with data. Usual 
parametrization of the PDFs at Qo reads 

xqi{x,Qo) = x^-'il - xf\ 

For this parametrization the behaviour of q at low x is motivated by the Regge phenomenol- 
ogy (see, e.g., book [^) and at high x, by the guark counting rules If such a simple 
form is insufficient for the fair data description, polynomial-like factors are added. Value of 
Qo is arbitrary, but it is natural to choose it as O(GeV) to allow for simple interpretation 
of the boundary PDFs. Meanwhile, it was recently shown in Ref. that the choice of 
Qo is important to provide the results stable with respect to the account of higher order 
QCD corrections (see also Ref. ||22|). At low Qq the twist 4 contribution extracted from the 
data is less sensitive to the choice of the renormalization scale /iR in Eg. ( p.6|) , than at high 
Qo- The as(Mz) behaves in opposite way, and then the choice Qq = 9 GeV^, made in our 
analysis, provides stability of the a^ and PDFs values if the NNLO QCD corrections are 
considered. 

Despite of the fact, that the evolution eguations have been used in the DIS data analysis 
for many years, no unigue approach for solving them exists. Analytical expressions can be 
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obtained only for the simplified splitting functions, and direct numerical approaches demand 
threefold integration, which is time consuming. There are semi-analytical approaches, based 
on expansion of PDFs in terms of selected sets of functions, but such approaches, as a rule, 
lead to loosing of the universality with respect to the choice of splitting functions and require 
careful control of the calculations precision. Due to the form of the evolution equation kernel 
is rather complicated, correct implementation of a sophisticated integration algorithm meets 
the difficulties. In the comparative analysis of different codes, used for the DGLAP equations 
integration, the codes of the CTEQ and MRST collaborations were found to contain the 
bugs (see Ref. |^). Taking into account these points, we use in the analysis our own code 
for direct numerical integration of Eq. (p.6| ), based on the Euler predictor-corrector algorithm 



(see Ref. (^g). This code allows one to modify kernels of the evolution equations in order 
to debug the code, to control the calculation precision, to take into account effects of new 
physics, and to implement special cases of evolution. Integration region can be expanded 
easily, and the integration precision is regulated by the external parameters of the code. 
For typical values of these parameters the code integration precision, as estimated using 
benchmark described in Ref. is given in Fig. |I|. One can see that the relative precision 
is better, than 0.001 in the region x ^ 0.5 and better, than 0.01 in the region x ^ 0.5. This 
is well enough for our purposes, since the errors on data are larger, than the integration 
errors for all x. 

Since Eq.( p.6|) is valid for massless partons only, the heavy quarks contribution, which 
is significant at low x, should be considered in a peculiar way. In the approach described in 
Ref. [^] the heavy quarks are considered as massless ones. They are included into the general 
evolution starting from a threshold value of Q, which is proportional to the quark mass, while 
at the values of Q lower the threshold these distributions are put to zero. Evidently, in this 
approach the heavy quarks contribution is overestimated in the vicinity of the threshold. 
Alternative way to consider the heavy quarks contribution is to calculate it using the photon- 
gluon fusion model of Ref. [^]. At high Q and low x "large logarithms" arise in the 
elementary cross section of this process, that may demand its resummation ||39|. Meanwhile, 
as it was shown in Ref. |40], the region of x and Q, where the resummation is really needed 
lays outside the region of the available DIS data. For this reason we calculated the c- and 
b- quarks contributions to the structure functions F2^l using the photon-gluon fusion model 
with the NLO coefficient functions of Ref. and the renormalization/factorization scales 



equal to y + 4m^ ,3 at the quark masses rric = 1.5 GeV and mb = 4.5 GeV. 

The LT contribution to the DIS structure functions is rather well understood both from 
theoretical and experimental points of view. Since this contribution depends on Q weakly 
one can reject the low Q data points and leave the data set, which is both statistically 
significant, and can be analysed within perturbative QCD in order to determine the LT x- 
dependence. The HT contribution is worse known, than the LT one. The theoretical analysis 
of the HT x-dependence is equally difficult as for the LT x-dependence and, as a result, it 
should be determined from data. Meanwhile, due to the fast fall of the HT contribution 
with Q it is sig nificant for < 10 GeV^ only. At very low Q the subtraction of the LT 
contribution, as calculated in perturbative QCD is problematic due to the rise of as. As a 
result, only the data for ^ 1 GeV^ Q can be used for the HT extraction and the results 
precision is poor. 
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Study of the possibility to separate the HT and LT contributions has a long history (see 



Refs. [^2H 45[]). Despite of that the Q-dependences of these contributions are different, in the 



limited range of Q the HT power corrections can simulate the logarithmic LT behaviour 



Moreover, as it was shown in Refs. , the power corrections can almost entirely describe 

the scaling violation observed for the DIS data, if the data precision is limited. In particular, 
this causes large correlations between the fitted values of as and the HT contribution. This 
correlation leads to the rise of the fitted parameters errors. The rise of errors is unpleasant 
effect, moreover, the fitted model non-linearity can become essential as a result. Besides, 
the fit results become less stable with respect to the change of the non-fitted parameters and 
adoptions of the fitted model, i.e., the theoretical errors on the fitted parameters rise also. 
Finally, if large correlations between parameters occur, the second derivative matrix for the 
minimized functional is poor determined and the calculations inaccuracies increase when its 
inversion. For this reason in order to get satisfactory precision of the parameters errors one is 
to guarantee better precision of the fitted model calculation, which may be time consuming, 
if manifold integration is involved. Due to this is the case for our analysis, estimation of the 
correlation coefficients between the fitted values of as and the HT contribution is non-trivial 
problem. 

III. DATA USED IN THE FIT AND STARTING PDFS 

We fit the PDFs to the data on the charged leptons DIS off proton and deuterium given 
in Refs. The data points with < 2.5 GeV^ were not used in the analysis in 

order to reject the region, where ag is rather large and the NNLO order QCD correction 
may be important. The points with x > 0.75, for which the nuclear corrections are large, 
were removed also. The data used in the analysis occupy the region 10~^ ^ x < 0.75, 
2.5 GeV^ < ~ 5000 GeV^. The number of data points for each experiment is given in 
Table |. 

The starting PDFs were initially parametrized at Qo = 9 GeV^ as follows: 

xq,{x,Qo)=Aix''^{l-xf'{l + YiV^ + lix) (3.1) 

and then the parameters 7, which agree with zero within errors, were by turn fixed at 
zero till such parameters existed. Evidently, the fit quality could not get worse, when such 
parameters are fixed. The PDFs functional form resulted from this simplification and used 
in the final fit reads: 

2 

XUy{x,Qo) = j^X^^il- xY-{l+-f^x), 



As 

xus{x,Qo) = —ri^x'''"'{l - x) 



xdy{x,Qo) = ^x"<^(l 
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xss{x,Qo) 



ATS 



r]sx"-^^{l - xf'% 



xG{x, Qo) = Agx'g (1 - (1 + + ^G^) 



where u, d, s, G are the up, down, strange quarks, and gluons distributions respectively; in- 
dices V and S correspond to the valence and sea quarks. The parameters , Nj and Aq 
were not fitted, instead they were calculated from the other parameters using the conserva- 
tion of the partons momentum and the fermion number. The parameter was calculated 
using the relation 



It is well known, that the charged leptons data do not allow to confine the sea quarks 
contribution. For this reason the parameter rjs was fixed at 0.42, which agrees with the 
recent results of the NuTeV collaboration, given in Ref. The other sea distributions 

parameters were constrained as Osu = ctsd = ctss, = (&su + &sd)/2. 

The DIS cross sections calculated from the QCD evolved PDFs using Eq.( p.l[ ) with the 
account of the TMC corrections given by Eq.( |2.3|) and the twist 4 contribution in additive 
form as in Eq.( p.2| ), were fitted to the cross section dataQ. The HT contributions to the 
proton and neutron structure functions F2 were parametrized by separate functions iff and 
H2 respectively, and the HT contributions to the proton and neutron structure functions 
Fl, by the common function H^, since the latter coincide within errors. The functions -fff '" 
and were parametrized in the model independent way: at x = 0., 0.1, . . . 0.8 their values 
were fitted to the data, and between these point the functions were linearly interpolated. 
The common approach for the PDFs global fits is to use data on F2 instead of the data 
on cross sections. Within this approach one ignores the fact, that the F2 values given by 
different experiments are often extracted from the cross sections using different values of Fl. 
In our fit the Fl contribution to the cross section was calculated iteratively and, efficiently, 
the data were reduced to the common value of Fl. Since the Fl contribution rises with y, 
the effect of this reduction is more important at high y. Due to the collision energy of each 
experiment is limited, the highest values of y correspond to the minimal values of x. For this 
reason the F2 data points shifts due to the reduction to the common value of Fl are not very 
significant in average, but at the edges of the experiments data regions may reach several 
percents. Note, that at low x the Fl value strongly depends on the gluon distribution and, 
hence, in the fit to the cross sections data an additional constraint for the gluon distribution 
occurs, i.e. it is better confined, as compared to the fit to the F2 data. 



^Since the high Q data from the HI and ZEUS experiments were corrected for the Z-boson 
contribution, Eg. (^H]) is applicable for these data also. 



2 / x[us{x, Qo) + ds{x, Qo) + Ss{x, Qo)]dx = As. 
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The TMC correction is most important for the SLAC data, less important for the BCDMS 
data, almost unimportant for the NMC data, and negligible for the others. Note, that our 
TMC correction to F2 given by Eq.( |2.3| ) differs from that applied in Ref. where the 
substitution 

Q) = F^^ii^uc. Q) (3.2) 

was used to account for the target mass effect. The numerical difference between these two 
approaches is maximal at high x and low Q, e.g., for the SLAC data it reaches 40%. Besides, 
our TMC correction, contrary to that given by Eg. (|3.2|), changes sign at x ~ 0.5. 



The deuterium data were corrected for the Fermi motion effect as in the model of Ref. 150 



with the deuterium wave function from Ref. [0]. The deuterium correction value rises with 



X and reaches 16% for the SLAC data. This correction was calculated iteratively in the fit 
to provide consistency of the analysis. The two-dimensional integrals involved in the model 
were calculated using the code of Ref. , which provides better numerical stability, than the 



standard codes based on the Gauss integration algorithm. For the calculations time saving 
we adopted, that the Fermi motion correction for the structure function xFi = F2 — Fl is the 
same, as for the structure function F2 (we checked that this adoption does not significantly 
affect the results). 



IV. FITTING PROCEDURE AND RESULTS 

The fitted parameters including the PDFs parameters, the value of as, and the coefficients 
of the functions if2,L were determined from the minimization of the functional 

= E (f^ - ^Ky^)EM - ex%), (4.1) 

K,i,j 

where Eij is inverse of the covariance matrix Cij, 

Cij = ^l^ijCTiaj + fifjifjf ■ r/f ), 

index K runs through the data subsets corresponding to the different experiments and 
the different targets within one experiment, indices i,j run through the data points in 
these subsets. The other notations: i/i are the measurements; cTj are the statistical errors; 
C,K are the renormalization factors; fi are the fitted model calculations depending on the 
fitted parameters; iff are the systematic errors vectors (the dimensions of these vectors 
for each experiment are given in Table | as NSE). The systematic errors were considered 
as multiplicative, that is natural way for the counting experiments. All systematic errors, 
excluding the normalization errors on the old SLAC experiments, were accounted for in the 
covariance matrix. The data from the old SLAC experiments, as they were given in Ref. 
|T[], are the result of re-analysis of the original experimental data published earlier (for the 



details see Ref. [Q). One of the purposes of this re-analysis was to renormalize the old 
data on the data from dedicated experiment SLAC-E-140. However, due to the latter did 
not release the proton target data, the renormalization of the proton data was performed 
using the experiment SLAC-E-49B as a "bridge" . Such technique certainly brings additional 
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uncertainties on the re-analysed data. In order to escape those uncertainties we performed 
the independent renormahzation of the old SLAG experiments without a "bridging" , that is 
possible in our case, since we use more proton data, than in the analysis of Ref. For this 
purpose we fitted the factors for each target of each old SLAG experiment independently. 
Alongside, to keep the analysis consistency the errors due to normalization uncertainties 
of the old SLAG experiments, given in Ref. 0], were cancelled out. For other experiments 
the parameters E,k were fixed at 1. The asymmetrical systematic errors on the ZEUS data 
were symmetrized, when including in the covariance matrix, and systematic errors on the 
BGDMS data for the proton and deuterium targets were considered as perfectly correlated. 
The statistical properties of the estimator based on covariance matrix (GME) were con- 



sidered in Ref. |TTI] in comparison with the statistical properties of the simplest estimator 



(SGE), based on the minimization of the functional 

2 _ V- (/i - Vi)' 



X, / y 9 1 

which is often used in particle physics for the analysis of data including the correlated ones 
as well. For the GME the fitted parameters systematic errors due to the data systematic 
errors are automatically included in the total error; for SGE the parameters systematic 
errors are estimated as the shift of the parameters under the shift of the data by the value 
of their systematic errors. The SGE dispersion is always larger, than the GME dispersion 



and, as it was shown in Ref. [0, the ratio of these dispersions can reach several units for 
realistic cases. It was shown also, that the GME is unbiased if the systematic errors on the 
parameters are not much more, than the statistical ones. In order to control the estimator 
bias one can trace the value of the net residual R, equal to the mean of weighted residual 
U ~y)l y^o"^ + {fvY ■ The values and the net residuals for the total data set and for each 
experiment separately calculated at the parameters values fitted using the GME are given 
in Table |. One can see, that the net residual value is within its standard deviation^ and 
the data description is good, excluding description of the ZEUS data. For more detailed 
analysis of the confidence of the ZEUS data description we calculated for those data the 
diagonalized residuals using the relation 



N 



where indices i, j run through the data points. If data are well described by fitted model, then 
for large N the values of rf obey the normal distribution, i.e. the Gauss distribution with 
zero mean and the dispersion equal to 1. The distribution of for the ZEUS data is given 
in Fig. 1^. Evidently it does not agree with the normal distribution, that is not surprising, 
since the data description is poor. Note, that the diagonalized residuals mean is small for 
the ZEUS data (0.05), meanwhile, the dispersion is equal to 2.1, i.e. it is far from the normal 
distribution dispersion. It is difficult to ascribe this discrepancy to the shortcoming of the 



The R standard deviation was calculated using Eq. (3.11) from Ref. |11|. 
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fitted model, since, as it seen from Fig. analogues distribution for the HI data agrees 
with the normal distribution perfectly, whereas both experiments gained similar statistical 
and kinematical coverage. One more possible explanation of this disagreement is, that the 
systematic errors given by the ZEUS collaboration are underestimated, but still are Gaussian 
distributed. In such cases the PDG scales the errors so that x^/NDP becomes equal to 1 



(see review [^). In our case this approach cannot be used, since number of independent 
sources of the systematic errors in the ZEUS experiment is large and a lot of variants of such 
rescaling can be applied. Besides, the distribution of residuals would remain non-Gaussian 
after the errors rescaling (see dashed curve in Fig. Driven by this consideration one can 
suppose that systematic errors on the ZEUS data are non-Gaussian distributed (but with 
zero mean) and then x^/NDP must not be equal to 1. If so, the fitted parameters, which are 
confined by the ZEUS data, also may be distributed in arbitrary way (see in this connection 
Ref. @])- -Due to exact estimation of the confidence intervals for unknown distribution is 



impossible, we recommend for this purpose, in particular for evaluating the PDFs errors at 
low X, the robust estimate of the confidence intervals, based on the Chebyshev inequality 
(see discussion in Ref. |rT| ). 

The dispersion of the net residual R is maximal for the SLAC-E-140, BCDMS, and 
NMC data sets. (Remind, that this dispersion rises with the increase of data correlation, 
full correlation corresponds to the dispersion of R equal to 1). Thus, one can conclude, 
that the account of the BCDMS and NMC data correlations has the largest impact on the 
analysis results, since number of points in the SLAC-E-140 data set is small. This conclusion 
is in line with the results of Ref. , where it was obtained, that in the combined fit to the 
non-singlet SLAC-BCDMS data the account of the BCDMS data correlations leads to much 
more significant shift of the parameters, than the account of the SLAG data correlations. 
The value of R for the total data set is well within its standard deviation, that confirms the 
fit unbiasness. 

The fitted PDFs parameters are given in Table 0. We underline, that in our fit the 
universality of the valence u- and d-quarks behaviour at low x is not initially assumed, 
contrary to the popular global fits practice, and the fit results confirm this universality with 
the few percents precision. At the same time the Regge phenomenology prediction (see, e.g., 



book |3T1) 

flu = Od = 0.5 (4.2) 

is in disagreement with the fit results^. A possible interpretation of this disagreement is, 
that Eq. ( |4.2|) , as it is deduced, is not related to a particular value of Q, while the QCD 
evolution does change the PDFs x-behaviour. As it was shown in Ref. ||5^, for the non- 



singlet distributions at low x this change is not very significant, but at least partially it 
can help to explain the observed disagreements. The values of the parameters Cu and 
agrees with the results of our earlier analysis of Ref. |^ and with the value of the parameter 
describing the low x-behaviour of the non-singlet neutrino structure function xF^, which 
was obtained from the fit to the CCFR data in Refs. lEBIBTll . For the obtained values of the 



^We especially underline this point, since Eq.(|4.2|) is often used for theoretical estimates. 
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parameters, which describe the valence u- and d-quarks behaviour at high x, the relation 
&d = Qu + 1 holds with good precision, in line with the quark counting rules. Meanwhile, 
the absolute values of these parameters deviate from the quark counting rules predictions 
6u = 3, 6d = 4. This disagreement can also be due to the QCD evolution, moreover for the 
non-singlet distribution the evolution effect is stronger at high x. 

As one can see from Table |I|, the systematic errors on the parameters describing the va- 
lence M-quark distributions at high x and the sea quarks distributions at low x are especially 
large. At the same time the ratio of the total error to the pure statistical one is 0(1) for any 
fitted parameter, that guarantees their unbiasness. In order to estimate the sensitivity of 
the parameters values to the approach used for the account of the systematic errors we per- 
formed the fit using the SCE and the fit with the statistical and systematic error combined 
in quadrature. Results of these fits are also given in Table |1[ One can see, that in the SCE 
fit the central values of some parameters are shifted by more, than two standard deviations, 
as compared with the CME fit and the shift is larger for the parameters with large ratio 
of the systematic errors to the statistical ones, in particular, for 6u and 7^. Nevertheless, 
the SCE fit provides correct estimate of the parameters, the only shortcoming of the SCE 
is that the SCE errors may be several times larger, that the CME errors. In our analysis 
maximal ratio of these errors is about 5 and within the errors the results of both fits agree. 
At the same time the fit with the statistical and systematic errors combined in quadrature 
does may give incorrect estimate of the parameters, since the data correlation information 
is lost in this case. As a result, the central values of some parameters, in particular, 6g and 
Osd, obtained from this fit are shifted from the CME fit results by the statistically significant 
values (see Table Some parameters errors obtained in these two fits are very different 
also, e.g., the errors on a^^Mi) and the parameters describing the gluon distribution at high 
X. These differences evidently may lead to the fake disagreements with another experimental 
results and cause discussions on new physics manifestation, if the results of the fit performed 
without the account of the data correlations are used for the comparison. (The example of 
resolution of such "disagreement" encountered in the comparison of the SLAC-BCDMS and 
LEP data on was given in Ref. [^ ). 

V. THE EXPERIMENTAL PDFS UNCERTAINTIES 

The fitted PDFs with their experimental errors, including both statistical and systematic 
ones are given in Fig. ^ and the relative experimental errors on the PDFs are given in Fig. ^. 
To estimate the separate contribution of the systematic errors to the total ones we calculated 
the parameters dispersions keeping the central values of the fitted parameters, but without 
the account of systematic errors on the data. Then we extracted these reduced dispersions 
from the total dispersions of the parameters and took the square roots of these differences 
as the systematic errors on the parameters. The ratio of systematic errors on the selected 
PDFs to their statistical errors is given in Fig. |. As it was noted above, the systematic 
errors impact is more important for the M-quark distribution at all x in question and for 
(i-quark distribution at low x. The PDFs errors, as well as their parameters errors, depend 
on the approach used for the accounting of systematic errors. The PDFs errors obtained in 
the CME and SCE fits are compared in Fig. ^ and one can see, that for the latter the errors 
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are several times larger generally. The errors on PDFs obtained from the CME fit in our 
earlier analysis of Ref. [0 are also given in the same figure. In that analysis we used data of 
Refs. with > 9 QeV^ and W > 4 GeV. At small and moderate x the errors on 

PDFs obtained in the earlier analysis are several times larger, that the PDFs errors obtained 
in the present analysis. At high x these errors are of the same order, and for some PDFs the 
earlier analysis errors are even smaller. This occurs due to in the analysis of Ref. the HT 
contribution was fixed at zero, that decreased the PDFs errors. The correlation coefficients 
matrix for the PDFs parameters is given in Table |T| and the selected PDFs correlation 
coefficients are given in Fig. [7|. The correlations are larger for the valence and sea quarks 
distributions. This can be readily understood, since these distributions contribute to the 
charged leptons DIS structure functions as the sum and hence can be separated hardly. Due 
to the large correlations between some PDFs the ratio of the systematic errors on their 
linear combinations to the statistical errors on these combinations may be not proportional 
to such ratios for the PDFs themselves. For example, as one can see in Fig. ^, the relative 
systematic errors on the sum of non-strange quarks distributions at low x are significantly 
smaller, than for the u- and (i-quarks distributions separately. 

The relative experimental error on the gluon distribution rises with x due to rapid falloff 
of the distribution itself. The prompt photon data were often used to better confine the gluon 
distribution at high x, but the prompt photon production data, which appeared recently, 
turned out to be in disagreement with the earlier data (see review [^). Besides, it was 
shown, that in the theoretical analysis of this process large uncertainties occur (see review 
||59||). For these reasons one cannot use the prompt photon data for pinning down the gluon 
distribution in a consistent way. In our analysis the gluon distribution at low x is determined 
by by the slope of the structure function F2 on Q (see Ref. and at high x, from the 
partons momentum conservation. The experimental errors on the sea quarks distributions 
are also rather large, since we did not use in the analysis the Drell-Yan process data. 

Unfortunately, the obtained PDFs and their errors suffers from definite model depen- 
dence. For example, if one releases the constraint Osu = dsd = Oss, the quarks distributions 
errors at low x rise significantly. Analogous effect is observed, if more polynomial factors 
are added to the starting PDFs. Such model dependence is inevitable, since it is impossible 
to determine a continuous distribution from limited number of measurements without addi- 
tional constraints. The model dependence is stronger for the PDFs correlated with another 
PDFs, e.g., for the sea and valence quarks distributions, while the model dependence is weak 
for the sum of these distributions. The gluon distribution is also insensitive to the variation 
of the quark distributions due to rather weak correlation with the latter (see Table |T|). 

The MRST and CTEQ PDFs are given in Fig. ^ in comparison with ours, although the 
comparison is incomplete, since the errors on the MRST and CTEQ PDFs are unknown. 
Note, that the difference between the MRST and CTEQ PDFs almost everywhere is smaller, 
than our PDFs errors. At high x it may occur due to those collaborations use in the analysis 
more data, but more probable explanation is that the MRST and CTEQ collaborations 
get similar results due to they use similar data sets. In particular, this means, that the 
difference between the MRST and CTEQ PDFs cannot be used as the estimate of the 
PDFs uncertainty. In the whole, with the account of our PDFs errors, there is no striking 
disagreement of our PDFs with the MRST and CTEQ ones. Our gluon distribution is slightly 
higher, than the MRST one at low x, but this disagreement is statistically insignificant. 
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Excess of our sea quarks distributions over the MRST and CTEQ ones at low x is statistically 
significant, but there are several reasons for it. Firstly, both collaborations use massless 
scheme for the account of the heavy quarks contribution, that can lead to the overestimation 
of this contribution, and the corresponding underestimation of the light quarks contribution 
at low X. Secondly, the MRST and CTEQ collaborations use in the analysis the CCFR 
neutrino data of Ref. |6l[], which confine the sea quarks contribution and which were recently 
corrected by the authors just at low x (see Ref. Finally, the discrepancy between the 

MRST and CTEQ PDFs is of the order of discrepancy between those PDFs and ours, i.e. 
one needs to perform a detailed analysis of all parametrizations to clarify this discrepancy. 
Excess of the u- and rf-quarks distributions over the MRST and CTEQ ones at x ^ 0.3 
is most statistically significant. We checked, that that this excess occurs due the MRST 
and CTEQ collaborations renormalize the BCDMS data by 1-2% downward. Since we do 
not apply such renormalization, our parametrization for F2, as well as the u- and (i-quarks 
distributions, lays higher. Besides, we applied the TMC correction and the correction on 
the Fermi-motion in deuterium, that also leads do the rise of the quarks distributions at 
moderate x. Note, that this excess may help to explain the excess of the TEVATRON jet 
production cross section data at transverse energies of = 200 — 400 GeV over the QCD 
predictions, since this cross section gets large contribution from the quark-quark scattering 
at X ~ 0.2. 



The comparison of our PDFs errors with the errors on PDFs of Ref. |jT^ is given in 



Fig. One can see that, despite of that in the analysis of Ref. [|T0[ an additional NMC data 
on the neutron and proton structure functions ratio and the CCFR neutrino data are used, 
our PDFs errors are smaller generally. We ascribe this difference to that in the analysis of 



Ref. |10] the SCE was used in the fit. This conclusion is supported by the comparison of 



the structure function F2 band, calculated from the PDFs of Ref. |TD[ , with the data used 
in that fit. The comparison is given in Fig. ^ One can see that the most left point error is 
smaller, than the error on the F2 parametrization of Ref. |TD[ for this point, i.e. SCE applied 



for that analysis uses information given by this measurement inefficiently. The qualitative 
explanation of such behaviour of the SCE is that for this estimator the fitted parameters 
systematic errors are basically determined by the data points with the largest systematic 
errors. The CME used in our analysis is more efficient, than SCE and, as one can conclude 
from Fig. |^, our error on the F2 parametrization is basically confined by the point with the 
lowest systematic error. The difference of the SCE and CME PDFs errors is more the more 
is the relative contribution of the systematic errors to the total one. As a consequence this 
difference is especially large for the w-quark distribution and it is demonstrative, that the 
error on the -u-quark distribution of Ref. |]TU| almost coincide with the w-quark distribution 
errors obtained from our SCE fit (see Fig. The error on d- and u-quarks distributions 
ratio at high x given by our PDFs is also smaller, as compared with this error given by the 
PDFs of Ref. |T0[ (see Fig. |). 

VI. THE THEORETICAL UNCERTAINTIES 

The theoretical uncertainties inherent for a phenomenological analysis cannot be ulti- 
mately defined, since in the study progress the set of such uncertainties may increase or 
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decrease. In our analysis we accounted for the following sources of the theoretical uncer- 
tainties: 



MC - the change of the c-quark mass by 0.25 GeV; 

SS - the change of the strange sea suppression factor by 0.1, in line with the estimate 
given by the NuTeV collaboration 



TS - the change of the heavy quarks threshold values from m^b to V 6.5mc,b; in accordance 
with the consideration of Sec.|l|; 

RS - the change of the renormalization scale in evolution the equations from Q/2 to 2(5; 

DC - the change of the deuterium nuclear model based on the account of Fermi-motion ||5^ 



on the phenomenological model of Ref. |6^. In view of the discussion of Refs. p^j65 



on the applicability of the model of Ref. to the light nuclei, one may suppose that 



this change leads to the overestimation of the corresponding error. 

These changes were made in turn and the fitted parameters shifts for each change were 
taken as the theoretical errors on the parameters. Sometimes in other similar analysis the 
PDFs theoretical errors due to the as and HT uncertainties are estimated using the same 
approach. In our analysis these errors are included in the total experimental errors, since 
both Q!s and the HT contribution are fitted. We underline, that the scales of the considered 
theoretical errors are rather conventional, since they are based on the "reasonable" estimates 
of the model uncertainties. For this reason the theoretical errors should be accounted for 
with certain cautions. 



VII. THE as VALUE AND THE HT CONTRIBUTION 

We obtained from the fit the value a^^Mz) = 0.1165 ± 0.0017(stat + syst). The exper- 
imental error on as obtained in our analysis is two times less, than in the NLO analysis 
of similar data set described in Ref. |2^, where the value as(Mz) = 0.1160 ± 0.0034(exp) 



was obtained. The contributions of separate sources of the theoretical errors on our value 
of as(Mz) are given in Table IV. One can see, that the largest contributions give uncertain- 
ties of the QCD renormalization scale and the heavy quarks threshold values (especially for 
6-quark). Combining all these contributions in quadrature, we obtain 



as(Mz) = 0.1165 ± 0.0017(stat + syst) i^;^^^^ (theor), (7.1) 
which agrees with the modern world average as(Mz) = 0.1184 ± 0.0031 given in Ref. []66| . 



Our estimate of the value is insensitive to the complication of the PDFs form, since it is 
almost uncorrelated with the PDFs parameters, in particular, with the gluon distribution 
ones (see Table fU^ . 



As it was recently reported in Ref. the net partons momentum for the PDFs, 



obtained from the data set similar to one used in our analysis, is not equal to 1, if one 
does not cut the data with ^ 10 GeV^. In particular, the net partons momentum 
obtained from the analysis of the world charged leptons DIS data with > 3 GeV^ is 
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< X >~ 1.08 ± 0.02, as it is given in Ref. pi]. The conclusion drawn from this observation 



is that the DIS data at low Q are irrelevant for the NLO QCD analysis and reliable results 
can be obtained from the fit to the data with > iq GeV^ > 10 GeV^ only. The 
value of as(Mz) = 0.114 ± 0.002 obtained in this analysis differs from ours. In order to 
perform comparison with this result, we repeated our fit without imposing the momentum 
conservation constraint on the PDFs and obtained that at = 9 GeV^ the net partons 
momentum is < a; >= 0.979 ± 0.029, which agrees with 1 and differs from the results of 



Ref. ||2^. For this reason we cannot support the conclusion of Ref. about irrelevance 



of the low Q charged leptons DIS data for the NLO QCD analysis. For more detailed 
comparison we performed the test fit with the cuts of Ref. and also obtained the lower 
value as(Mz) = 0.1098 ± 0.0055, but with the error, which is significantly larger, than one 
obtained in Ref. [Q, and the as value obtained in this test fit is in agreement with ( [7.1| ) 
within the errors. The observed difference of the errors evidently occurs due to in our 
analysis we simultaneously fit both the ctg value and the HT contribution to F2. As it was 
shown in Refs. p6|JT6[l , the latter are strongly correlated, that certainly leads to the rise of 
the parameters errors. In support of this conclusion, if in our test fit the HT contribution 
is fixed, the as error falls from 0.0055 to 0.0014. However, the results of the fit with the 
HT fixed are model dependent and essentially the decrease of the experimental error is 
accompanied by the uncontrolled rise of the theoretical errors. 

The HT contributions to the nucleon structure functions Fl and to the proton and 
neutron structure functions F2 are given in Fig. |T0| and in Table 0. It is interesting that up 
to minimal x the twist 4 contribution to the structure function F2 is non-zero, that coincides 
with the results of Ref. on the analysis of the NMC data. The deviation of the Fl twist 4 
contribution off zero at low x is even more significant. As one can see from Table 0, the HT 
contributions to F2 and to Fl at low x are very sensitive to the approach used to account for 
the systematic errors on data. This is due to at low x the HT contributions are determined 
from the comparison of the data at the kinematical edges of different experiments, where 
the systematic errors are largest as a rule. Note, that the HT parameters errors obtained 
in the CME fit are 2-3 times smaller, than in the SCE fit, as well as the PDFs parameters 
errors. 

The twist 4 contributions obtained at the different values of the QCD evolution equations 
renormalization scale /xr are given in Fig. ^ The evident dependence of H2 on /xr at low 
X indicates that in this x-region the twist 4 contribution to F2 can simulate the effect of the 
NNLO corrections to the splitting functions P. Analogous effect for the structure function 
xF^ was demonstrated in Ref. [^] , while the direct observation of the re-tuning of the twist 
4 contribution to xF^ due to the account of the NNLO corrections was reported in Ref. [f^l| . 
At the same time the /xr dependence of i^L and of H2 at high x is not so strong. The 



explanation of such behaviour is given in Ref. |3J]. As it was also shown there, due to 
the HT contribution can partially absorb the NNLO corrections effects, the /xr dependence 
of the as value obtained in the simultaneous fit of the HT contribution and as is weaker, 
than in the fit with the HT contribution fixed. In particular, due to this absorption, the 
tts renormalization scale error obtained in our analysis is smaller, than in the analysis of 
Ref. [gl. 

The difference of the HT contributions to the proton and neutron structure functions F2 
is given in Fig. [T^. One can see, that at low x these contributions coincide within errors. This 
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is in disagreement with the results of Ref. |^.In that paper the data on the difference of the 
proton and neutron structure functions F2 are compared with the calculations based on the 
standard PDFs and found to be lower than that calculations at a: ~ 0.3. This discrepancy 
was attributed to the existence of the large HT contribution to the difference of the proton 
and neutron structure functions F2. We do observe the statistically significant deviation of 
H2 — H2 off zero, but at a; ~ 0.7 instead of x ~ 0.3. Unfortunately, this difference strongly 
depends on the deuterium nuclear corrections model at large x (see Fig. |TT]) and in order 
to obtain a reliable estimate of the twist 4 contribution to F2 an additional comparative 
analysis of the deuterium models is needed. 



VIII. THE PARTON LUMINOSITIES AT THE FNAL AND LHC COLLIDERS 

All errors on the hard processes cross sections due to the PDFs uncertainties are con- 
centrated in the parton luminosities (PLs), defined as 



L,,(M) = - / —q,{x,M^)q,{T/x,M^ 

S JT X 



where s is the s.c.m. energy squared; M is the produced mass; r = M"^ / s\ i and j mark the 
parton species. Since the PLs errors strongly depend on the latter, one is to estimate the 
impact of the PDFs errors on the calculated cross sections errors in each particular case. 
Our PDFs total errors, comprised of the theoretical errors combined with the experimental 
ones are given in Fig. ^. Despite of that the data set used for the extraction of our PDFs 
is limited by the DIS data, the PDFs errors are rather small at low x, i.e. in the region 
especially important for the FNAL and LHC experiments. The valence quarks distributions 
errors are small at high x also (see Fig. ^). The experimental errors dominate for the sea and 
gluon distributions at high x only (note, that this is not the case for the SCE fit, as one can 
see from the comparison of Fig. ^ and Fig. As one can see from Fig. 0, the dominating 
source for the gluon distribution at low x is the RS uncertainty, for the sea distribution at 
low X, the MC one, for the (i-quark distribution, the DC one. Remind, that in our analysis 
the errors due to the uncertainties of the ag value and the HT contribution are included into 
the experimental error. To estimate their contribution to the total error we re-calculated 
the PDFs dispersions fixing the as value and the HT contribution by turn, then extracted 
obtained dispersions from the nominal dispersions calculated with these parameters released. 
The square roots of these differences were taken as the PDFs errors due to the as and the 
HT uncertainties respectively. The ratios of these errors to the total PDFs errors are also 



given in Fig. |T2[ One can see, that the ttg uncertainty affects the gluon distribution only, 
while the HT uncertainty contributes to the errors of all PDFs. 

The errors on the PLs relevant for the most common processes at the energy of the FNAL 
collider are given in Fig. 0. The upper limit of the pictures was chosen so that the PLs at 
the upper limit is ~ 0.01 1/pb, i.e. corresponds to the maximal sensitivity of the planned 
experiments. One can see, that at the FNAL collider energy the theoretical errors dominate 
over the experimental ones at M < 0.2 TeV and vice versa at M ^ 0.2 TeV. The total 
PLs errors for the FNAL collider generally do not exceed 10% at M ^ 0.2 TeV, while for 
the quark-antiquark PL the total error is smaller, than 10% almost for all M in question. 



17 



The PLs pictures for the LHC energy, given in Fig. |T^, approximately reproduce the FNAL 
pictures with the produced mass M scaled in 5 times and the quark- ant iquark PL replaced 
by the quark-quark PL. 

Due to the PDFs correlations generally are not small (see Fig. the account of these 
correlations may affect the calculated hard processes cross sections errors. In some cases the 
PLs errors may cancel in their ratio, as in the example given in Table 0. Calculating the 
theoretical errors on the hard processes cross sections one is also to take into account the 
correlations of PDFs with the elementary processes cross sections, if the latter depend on 
the parameters responsible for the PDFs theoretical uncertainties. Besides, the RS PDFs 
uncertainty may be compensated by the NNLO corrections to the elementary processes cross 
sections. 



IX. CONCLUSION 

Significant part of the studies planned for the next generation hadron-hadron and lepton- 
hadron colliders is devoted to the precise Standard Model checks (see, e.g., review |]TD|). Such 
studies certainly imply careful control of all possible uncertainties, including the PDFs er- 
rors. The PDFs obtained in our analysis are supplied by the experimental and theoretical 
errors and can be used for the correct estimate of the calculated hard processes cross sections 
uncertainties, necessary for a precise phenomenological comparison aiming to detect a man- 
ifestation of new physics (e.g., compositness in proton-proton and electron-proton collisions, 
the partons recombination at low x, precise determination of the W and Z masses, etc.). 
A particular feature of our PDFs is that they were obtained using efficient estimator and, 
as a result, have minimal errors. The convenient code allowing to account for the PDFs 
uncertainties in the Monte Carlo calculations is accessible through the computer network^. 
Using the current version of this code one can obtain the random Gaussian smeared PDFs 
values with the account of the experimental and the theoretical uncertainties and their cor- 
relations. The special parameters allow one to scale the dispersions corresponding to the 
separate sources of the PDFs uncertainties to give user the possibility to study effects of 
each uncertainty and vary the confidence level of the errors on the calculations results. 
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FIGURES 




FIG. 1. A is relative precision of our code used for the evolution equations integration. Indices 
u and d correspond to the valence quarks; S, to the sea quarks; G, to gluon. 





FIG. 2. The distribution of diagonalized residuals for the ZEUS and HI data (full curves: 
normal distribution, dashes: the Gauss distribution with the dispersion and the mean equal to the 
dispersion and the mean of the residuals distribution) . All curves are normalized to the number of 
points in each experiment. 



22 




FIG. 3. The la experimental error bands for our PDFs at = 9 GeV^ (full lines). For 
comparison the nominal MRST99 (dots) and CTEQ5 (dashes) PDFs are also given. 
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FIG. 4. The relative experimental PDFs errors [%] (full lines: the total errors, dashes: the 
experimental ones). For comparison the relative experimental errors on the PDFs of Ref. [^] are 
also given (dots). 




10 10 10 10 

FIG. 5. The ratio of the systematic errors on the fitted PDFs to the statistical ones (full lines: 
the gluon distribution; dashes: the total sea one; dots: the d-quark one; dotted-dashes: the li-quark 
one). 
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FIG. 6. The relative experimental PDFs errors [%] (full lines: our analysis, dashes: the analysis 
of Ref. [|lO|). For comparison the relative experimental PDF errors obtained in our analysis from 
the SCE fit are also given (dots). 
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FIG. 8. The la experimental error bands for F2 calculated using different PDFs (full line: our 
PDFs; dashes: the PDFs of Ref . Q). Circles: the SLAC data; squares: the BCDMS ones. 
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FIG. 10. The twist 4 contribution to the proton structure function F2 and to the nucleon 
structure function Fl (full circles: = Q; open circles: = 2Q; squares: //r = Q/2). For 
better view the points corresponding to different /iR are shifted to left-right along the x-axis. 




FIG. 11. The difference of the twist 4 contributions to the neutron and proton structure func- 
tions F2 obtained in the fits using the different deuterium models (full circles: the Fermi- motion 
model; open circles: the model of Ref. ||6^. For better view the points corresponding to different 
models are shifted to left-right along the x-axis. 
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FIG. 12. The ratios of the PDFs errors due to separate sources to the total PDFs errors (fuh 
lines: the gluon distribution; dashes: the non-strange sea one; dots: the d-quark one; dotted-dashes: 
the u-quark one). THEOR means the sum of the MC,SS,RS,TS, and DC contributions. 
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FIG. 13. The relative errors on selected PLs for the FNAL collider (full lines: experimental 

errors, short dashes: RS; dotted-dashes: TS; sparse dots: DC; dense dots: MC; long dashes: SS). 

Other notations: Lqq = L^^ -|- L^a + -^du; -^qq = -^ud + -^da; -^(q+q)G = -^uG + -^uG + -^dG + -^dG- 
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TABLES 



TABLE I. Total number of points (NDP), number of independent sources of systematic errors 
(NSE), x^/NDP and the net residual R for each experiment (standard deviation of R is given in 
parenthesis). Also the renormalization factors ^ for the old SLAC experiments are given. 





NDP/(1 










experiment 


proton 


deuterium 


NSE 


yVndp 

A. / ' -^-^ 


R 


SLAC-E-49A 


58 / 1.8 ± 1.3 


58 / -0.4 ± 1.2 


3 


0.52 


-0.05(0.23) 


SLAC-E-49B 


144 / 2.0 ± 1.3 


135 / -0 1 ± 1 3 


3 


1.20 


22f0 29) 


SLAC-E-87 


90 / 2.0 ± 1.2 


90 / 0.2 ± 1.2 


3 


0.91 


0.01(0.37) 


SLAC-E-89A 


66 / 4.2 ± 1.8 


59 / 1.2 ± 1.9 


3 


1.34 


-0.18(0.45) 


SLAC-E-89B 


79 / 1.5 ±1.2 


62 / -0.7 ±1.2 


3 


0.82 


0.46(0.49) 


SLAC-E-139 




16 / 1.0 ± 1.2 


3 


0.64 


-0.10(0.43) 


SLAC-E-140 
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4 


0.89 


0.51(0.86) 


BCDMS 


351 


254 


9 


1.15 


0.07(0.68) 


NMC 


245 


245 


13 


1.32 


0.05(0.62) 


HI (94) 




147 


5 


0.96 


0.11(0.25) 


ZEUS (94) 




188 


20 


2.14 


0.32(0.34) 


FNAL-E-665 


47 


47 


10 


1.23 


0.38(0.38) 


Total 


1080 


1327 


79 


1.20 


0.12(0.22) 
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TABLE II. The fitted Os and the PDFs parameters values (I: the CME fit; II: the SCE fit; 
III: the fit with the statistical and systematic errors combined in quadrature). Given errors on 
the parameters include both statistical and systematic errors, pure statistical errors are given in 
parenthesis. 







I 


II 


III 


Valence quarks: 














0.693 ± 0.033(0.027) 


0.715 ± 0.114(0.029) 


0.703 ± 0.035 






3.945 zb 0.050(0.039) 


4.119 ± 0.257(0.038) 


4.037 ± 0.049 




72 


1.29 lb 0.44(0.37) 


1.39 ± 1.86(0.40) 


1.42 ± 0.49 




ad 


0.725 ± 0.086(0.082) 


0.703 ± 0.172(0.094) 


0.717 ±0.13 




bd 


4.93 lb 0.13(0.12) 


4.83 ± 0.27(0.17) 


5.00 ± 0.17 


Gluon: 












ac 


-0.225 ±0.035(0.031) 


-0.169 ±0.065(0.029) 


-0.135 ± 0.044 




bo 


6.1 ±2.1(1.8) 


4.9 ±5.6(1.7) 


4.07 ± 1.3 




1? 


-2.63 ± 0.83(0.71) 


-3.41 ± 0.99(0.45) 


-4.06 ± 0.48 




1§ 


4.7 ±2.9(2.4) 


4.44 ±3.4(1.3) 


5.41 ±1.2 


Sea quarks: 












As 


0.166 ±0.011(0.0095) 


0.167 ±0.025(0.011) 


0.167 ±0.017 




«sd 


-0.1987 ±0.0067(0.0050) 


-0.1853 ± 0.0181(0.0050) 


-0.1833 ± 0.0075 




bsd 


5.1 ± 1.4(1.3) 


5.4 ±2.8(1.4) 


4.9 ±2.1 




Vn 


1.13 ±0.11(0.087) 


1.10 ±0.23(0.086) 


1.16 ±0.16 




bsu 


10.29 ± 0.97(0.81) 


10.56 ± 3.2(0.83) 


11.2 ± 1.1 




as(Mz) 


0.1165 ± 0.0017(0.0014) 


0.1138 ± 0.0044(0.0021) 


0.1190 ± 0.0036 
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TABLE III. The correlation coefficients for the starting PDFs parameters. The largest coeffi- 
cients are printed in bold. 





Ou 


feu 


72" 




fed 




Osd 


fesd 




fesu 


a-G fee 


1? 


72^ 




1.00 


























feu 


-0.84 


1.00 
























72" 


-0.97 


0.92 


1.00 






















Od 


-0.09 


-0.09 


0.05 


1.00 




















bd 


-0.21 


0.02 


0.19 


0.71 


1.00 


















As 


-0.14 


0.34 


0.24 


-0.86 


-0.54 


1.00 


















0.58 


-0.45 


-0.55 


0.37 


0.16 


-0.46 


1.00 














bsd 


-0.05 


-0.10 


0.00 


0.97 


0.54 


-0.88 


0.40 


1.00 














0.25 


-0.13 


-0.23 


-0.69 


0.24 


0.47 


0.01 


-0.78 


1.00 










bsn 


0.83 


-0.74 


0.86 


-0.14 


-0.16 


-0.24 


0.62 


-0.10 


0.44 


1.00 








«G 


0.23 


-0.22 


-0.23 


0.37 


0.20 


-0.38 


0.53 


0.37 


-0.21 


0.18 


1.00 






bo 


0.18 


-0.17 


-0.20 


0.11 


0.17 


-0.08 


-0.10 


0.06 


-0.11 -0.02 


0.27 1.00 






7f 


-0.36 


0.34 


0.36 


-0.45 


-0.30 


0.48 


-0.52 -0.44 


0.18 


-0.30 


-0.82 -0.47 


1.00 




72^ 


0.34 


-0.34 


-0.36 


0.28 


0.26 


-0.32 


0.15 


0.23 


-0.11 


0.20 


0.46 0.89 


-0.77 


1.00 


as(Mz) 


0.22 


-0.31 


-0.18 


0.01 


-0.05 


-0.05 


0.04 


-0.01 


0.04 


0.17 


0.01 -0.39 


0.03 


-0.18 



as(Mz 



1.00 



TABLE IV. The as(Mz) theoretical errors due to different sources. 

Source Value 

MC ±0.0003 

SS ±0.0001 

pq _l_0.0026 

="=0.0024 

TS -0.0020 

DC -0.0012 
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TABLE V. The fitted twist 4 contributions (I: the CME fit; II: the SCE fit; III: the fit with the 
statistical and systematic errors combined in quadrature). Given errors on the parameters include 
both statistical and systematic errors, pure statistical errors are given in parenthesis. 



X I II III 



0. 


-0.085 ± 0.026(0.020) 


-0.124 ±0.051(0.020) 


-0.132 ±0.035 


0.1 


-0.057 ±0.019(0.014) 


-0.107 ±0.076(0.014) 


-0.094 ±0.021 


0.2 


-0.024 ± 0.012(0.0097) 


-0.057 ±0.049(0.010) 


-0.054 ±0.016 


0.3 


-0.010 ± 0.010(0.0089) 


-0.027 ± 0.024(0.0090) 


-0.017 ±0.015 


0.4 


0.002 ± 0.010(0.0089) 


0.002 ± 0.024(0.0090) 


-0.002 ± 0.016 


0.5 


0.0292 ± 0.0085(0.0074) 


0.041 ±0.020(0.0079) 


0.025 ±0.015 


0.6 


0.0522 ± 0.0078(0.0069) 


0.068 ±0.017(0.0074) 


0.051 ±0.013 


0.7 


0.0535 ± 0.0061(0.0055) 


0.074 ± 0.013(0.0058) 


0.056 ±0.010 


O.o 


0.04o5 ± 0.0UD4(U.U0Dlj 


0.0545 ± U.UOoo^O.OUbUj 


0.0471 ± O.OUoo 


0. 


0.332 ± 0.046(0.033) 


0.13 ±0.11(0.033) 


0.028 ± 0.061 


0.1 


0.108 ±0.020(0.016) 


0.117 ±0.065(0.016) 


0.118 ±0.022 


0.2 


0.094 ±0.019(0.015) 


0.145 ±0.047(0.015) 


0.097 ±0.021 


0.3 


0.096 ±0.018(0.016) 


0.133 ±0.031(0.016) 


0.115 ±0.021 


0.4 


0.014 ± 0.017(0.015) 


0.040 ± 0.027(0.015) 


0.033 ± 0.019 


0.5 


0.0179 ± 0.0088(0.0068) 


0.023 ±0.014(0.0069) 


0.015 ±0.011 


0.6 


0.0031 ± 0.0094(0.0076) 


-0.016 ±0.024(0.0076) 


-0.0033 ± 0.0089 


0.7 


0.0195 ± 0.0064(0.0056) 


0.008 ± 0.016(0.0055) 


0.0134 ± 0.0067 


0.8 


0.024 ±0.012(0.012) 


0.01 ±0.023(0.012) 


0.012 ±0.014 


0. 


0.054 ±0.050(0.041) 


0.045 ±0.112(0.041) 


0.095 ± 0.077 


0.1 


0.031 ±0.027(0.026) 


0.041 ± 0.047(0.026) 


0.003 ± 0.037 


0.2 


-0.017 ±0.018(0.017) 


0.024 ±0.046(0.017) 


-0.014 ± 0.024 


0.3 


0.010 ±0.017(0.016) 


0.052 ± 0.038(0.016) 


0.014 ±0.021 


0.4 


0.023 ± 0.016(0.015) 


0.047 ± 0.037(0.015) 


0.036 ± 0.019 


0.5 


0.0068 ±0.011(0.010) 


0.009 ± 0.026(0.011) 


0.019 ±0.016 


0.6 


-0.029 ±0.0091(0.0086) 


-0.037 ±0.016(0.0092) 


-0.022 ±0.015 


0.7 


-0.052 ± 0.0073(0.0068) 


-0.073 ±0.014(0.0071) 


-0.055 ±0.011 


0.8 


-0.075 ±0.011(0.010) 


-0.079 ± 0.014(0.010) 


-0.073 ± 0.013 
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TABLE VI. The relative errors on the PLs involved in the calculations of the W and 
Z production cross sections at the FNAL collider {L^^ = L^^ + L^^, Lz = L^^ + -^dd' 

Lw/Z = (^ud + ^du)/ (-^uu + -^dd))- 



stat+syst RS TS SS MC DC 

ALwi%) 1.5 - - 1.2 1.1 1.5 

ALz(%) 1.2 - - 1.2 1.1 1.5 

^Lw/z{%) 0.7 _____ 



34 



